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1 Introduction 

As it is known, nonlinear systems are used in a variety of applications and have 
been the focus of research for a number of years. Many of mathematical tools have 
been used for the analysis and design of such systems. Among these tools, algebraic 
and computational methods have had a major effect. These tools have led to the 
development of new algorithms and techniques, which in turn have allowed design 
methods to be accomplished with greater speed and efficiency, [12]. 
Such computational methods have been applied in the study of modeling, the prob- 
lem of feedback linearization, and in the global optimization problem. ([1],[2],[12], 
to mention but a few). All of these problems were investigated mostly in relation to 
continuous systems, and, in some cases, in relation to non-linear discrete time sys- 
tems, [13]. By definition, a discrete time system evaluates input and output signals 
over a countable number of time instants. In some cases, a discrete time system is 
obtained from continuous time systems through sampling at certain time instants 
[15]. A further case concerns systems that are naturally and directly described in 
discrete form, typically in financial or economic systems, [7]. There is a rich liter- 
ature devoted to the study of discrete systems. Certain works approach the issue 
through analytical tools, [4], [5] and others by using algebraic methodologies, like 
differential algebra or rings theory, [1],[2], [3], [6], [23]. 

Factorization is a popular algebraic method dealing with linear as well as nonlinear, 
continuous or discrete, systems. In linear systems the central idea is that of factoriz- 
ing the transfer matrix as the ratio of two stable rational matrices, [15]. In nonlinear 
systems construction of coprime factorizations have been studied by several authors. 
This setting was mainly applied to systems that can be described by input-output 
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representations, [2] , [21] . 

The above factorizations can be used effectively in many design problems, [21], [24]. 
The simplification-linearization method is among them. This method can permit us 
to determine systems with a " less " complex structure, which are equivalent to the 
original nonlinear system, i.e. they will give the same output under the same input 
and identical initial conditions. Among these " simpler " systems the linear systems 
are the most desirable ones, because their behavior is well known. Therefore, we 
may ask when a nonlinear system is equivalent to a linear one. This is a version of 
the linearization problem. 

The present paper focuses on the simplification-linearization problem of nonlinear 
discrete input-output systems of the form: 

y{t) + J2 (^ivit + (^iii2y{t - 'i'i)y{t - ^2) h \- 

+ 1212 ■■ -12 (^hi2...iny{t - ^1) • • • y{t - in) = 

= J2 bjU{t - j) + --- + J2J2---J2 ^hh-jm'^it - ji) ■■■U{t- + 

+EE^«^(^-%(^-0 + --- + EEE---EE---E 

Ckik2...K,hi2...ivy{t - ■■■y{t- K')u{t - /i) ■ ■ ■ u{t - Q (1) 

Equation (1) transforms causal input signals (i.e. u{t) = for t < 0) to causal 
output signals. A set of initial conditions y(0) = yo, y{l) = yi, . . ., y{k) = yk is 
always assigned to (1) and the lowest delay output term (that is y{t)) appears in 
the linear part of the system. The systems of the form (1), which contain products 
among inputs and output signals that are sometimes called 'cross-products', encom- 
pass a broad variety of nonlinear discrete systems. We obtain these either through 
transformations of nonlinear discrete state-space representations into input - output 
forms [19], or when we use Taylor's expansion method to approximate other more 
general nonlinear discrete systems [14], [19]. These are employed in signal processing 
theory, whenever it is necessary to construct nonlinear representations of discrete 
signals, (they are an extension of the infinite impulse response filters to a nonlinear 
set up, [14]), in nonlinear time-series analysis and in adaptive control, in the context 
of designing nonlinear adaptive controllers, [22]. 

In order to describe the above systems algebraically we create an appropriate frame- 
work by using the so-called 5e-operators to deal with cross-products. This operator, 
introduced in [9], is an extension of the simple 5 and e-operators that are used in 
the algebraic description of nonlinear discrete input - output systems without cross- 
products [10], and which are also an extension of the simple shift operator q in linear 
discrete systems, [20]. Via these operators we can define the so-called 5,e and 5e - 
polynomials. By means of these polynomials we can formally rewrite (1) as follows: 
Ay{t) = Bu{t) + C[y{t),u{t)], where A a (^-polynomial, B an e-polynomial and C 
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a 5e-polynomial. Among these polynomials, we can define two product operations: 
the dot-product (denoted by '•'), which corresponds to the usual product among 
polynomials, and the star-product (denoted by '*'), which corresponds to the sub- 
stitution of one polynomial by another, or, in more system-oriented terminology, 
to the cascade connection of systems. Though there is a similarity between this 
algebraic background and others, there is also the interesting peculiarity that the 
set of 5e-polynomials with respect to the star-product is not a ring [1]. 
Two systems Aiy{t) = Biu{t) + Ci[y{t),u{t)] and A2y*{t) = B2v{t) + C2[y*{t),v{t)] 
are equivalent, if y{t) = y*(t), whenever u{t) — v{t) and initial conditions are 
identical. A " simplification " method for a given non-linear system consists in 
discovering systems that are equivalent to the original one but have a " less " complex 
structure. This procedure allows us to replace, if necessary, the original system with 
the simpler one. Clearly, of all such 'simpler' systems, linear systems are the most 
desirable ones. Accordingly, in this paper we develop a method for discovering 
linear systems which are equivalent to a given non-linear system. In other words, 
we are looking for a linear system Aiy*{t) = Biv{t), such that y*{t) = y{t), whenever 
u{t) = v{t). In other words we have a kind of the linearization procedure devoted 
to open-loop systems, thus, in our approach we do not use any feedback law design. 
To face the simplification-linearization problem we separate the systems of the form 
(1) into two categories, the first consists from nonlinear systems without cross- 
products, the second from nonlinear systems with cross-products. To work with 
systems of the first class we follow the next steps: First, be means of a symbolic 
algorithm, we factorize the nonhnear system as ^4 * Ly{t) — B* Mu{t), where L, M 
are linear 5 and e-polynomials correspondingly and A, B nonlinear. The novelty 
is that the coefficients of all the above polynomials are not constant number but 
functions of certain undetermined parameters Wij, Sij. By giving to these parameters 
certain values, we endow the polynomials with specific properties. Therefore, using 
certain methods of computational algebra we can find values of Wij, Sij so that A — B 
and hence the nonlinear polynomial can be simplified to the linear one Ly(t) — 
Mu{t). We have to mention here that this simplification procedure is something 
which is well known in the linear case. Indeed, let q be the simple shift operator, 
i.e. q~^y{t) = y{t — 1) and Ay{t) = Bu{t) a linear system. A, B polynomials of 
the variable q~^. If we re- write it in the form C ■ A'y{t) — C ■ B'u{t) and C 
is a stable polynomial, then, by performing a division, we take the linear system 
A'y{t) = B'u{t), which gives the same output with the original system, under the 
same input and initial conditions [20]. To deal with polynomials of the second 
class, that is those with cross-products, we work in an analogous way: First we 
factorize the system as A* [L, M] = 0, where L, M linear 5 and e-polynomials and 
A a (5e-polynomial with parametrical coefficients. Proper choice of values for the 
parameters may give to A a certain structure, the so-called homogeneous structure. 
If this happens, then we can write the nonlinear system as A * {Ly{t) + Mu{t)), A a 
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nonlinear polynomial, and then by simplification of A we take the linear equivalent 

system Ly{t) = —Mu{t). 

The particularities of our methodology are: 1) Its symbohc algorithmic orientation 
that allows a direct computer implementation. Actually, we create a proper algebraic 
framework so that specific algorithms to be applied. These algorithms permit us 
to transform questions from system theory, such as linearization, to computational 
algebra questions. Indeed, in order to take solutions to our simplification procedure 
we have to solve systems of polynomial equations, this can be done by means of 
certain algorithms of computational algebra. We note here, that our approach is 
symbolical and not numerical. 2) It gives not a single solution but a class of them. 
Indeed, if the above mentioned systems of polynomial equations do not accept a 
single solution but many of them, then to its such solution corresponds a linear 
system, equivalent to the nonlinear one. 

Finally, we have to make clear, that in the present paper we are not involved with 
internal stability questions. We are mainly interest to describe the simplification 
procedures formally. The text is divided into two parts. In Part I wc present the 
algebraic framework and the algorithms, in Part II we are dealing separately, with 
the simplification-linearization problem of systems without cross-products and of 
systems which include cross-products. Throughout the paper, N, Z, R will denote 
the set of natural, rational and real numbers, respectively. 

2 Part I - The Algebraic Preliminaries 

In this section we present the algebraic preliminaries, needed for the development 
of the simplification algorithms. 

2.1 The Multiindices 

Let A; be a positive integer. A subset of the set U^^j^Z" is called a set of multiindices 
and it is denoted by I. A set of multiindices may be finite or not. We denote the 
elements of I by i = {11,12, - ■ ■ ,in)- Usually, we put the elements of a multiindex i in 
an ascending way, that is ii < ^2 ^ " ■ ■ ^ "^n- The quantity ii + 12 + ■ ■ ■ + in , called 
the degree of i, is denoted by deg{\). Wc can define the following operations among 
multiindices. These will be used in the results in systems simplification below. 

• Let i = (ii, ^2, • • • , "^n) be an multiindex. An multiindex of the form j = 

), with G {ii, ^2, • • • , in} and cxi 7^ 7^ ■ ■ ■ 7^ cxfc, is called a 
A;-subindex of i. 

• Let i = {ii,i2, ■ ■ ■ ,in) be an multindex and j = (i<^i, i<^2! • • • ) Vfc) ^ subindex 
of i. Their difference is defined as follows: i j = {ixi,ix2, ■ ■ ■ ■, ixm)^ '^i^^ 

£ {^1) ^2) • • • ) in} {Vi' '^'P2i ■ ■ ■ 1 i<Pk} ■ 
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• Given two multiindices i = (ii, • • • , ^fc) and j = {ji, j2, ■ ■ ■ , jx), the new 
multiindex i©j is defined just juxtaposing j after i. Explicitly, i©j = 

Ji,«2,«3, J2, ■ ■ ■ ,«fc, Ja), where ii < ji < 12 < is < j2 < ■ ■ ■ < ik < j\- 

• Let j = (ji, J25 • • • ,jm) be a multiindex. We define the pointwise sum j+(p, (p 
a real number, as follows: j+ip = {ji + (p, j2 + <f, . . . ,jm + p)- 

Finally, a set I of multiindices can be ordered in a lexicographical way as follows: Let 
i = (^ij ^2j • • • , in) and j = (ji, j2, . . . , jm) be two multiindices, with elements ordered 
in an increasing way. We can say that the multiindex i is " less" than the multiindex 
j, and we write i < j, if either n < m 01 n — m and the right- most nonzero entry 
of the vector j — i is positive. Since for any two multiindices we can immediately 
deduce which multiindex is "larger" and which "smaller", the above lexicographical 
order is well defined. 

2.2 The (5e-operators 

We introduce the 5e-operator to deal with cross-products among input and output 
signals, which acts upon a pair of sequences. This initially appeared in [9] and 
has been examined in [17], [16]. Let y{t),u{t) be real sequences, defined over the 
set of integers Z and let F be the set of causal sequences. Thus if y{t),u(t) G F 
then y{t) = u{t) = for t < 0. Let i = (ii, ^2, • • • , ^m), j = (ii, ^2, • • • , in) be 
multiindices. The operator Si x 6^ : F x F ^-^ F is defined as: 6i x 6j[y{t),u{t)] = 
y{t—ii)y{t—i2) ■ ■ ■ y{t—im)u{t—ji)u{t—j2) • • ■u{t—jn)- This means that the operator 
5i acts exclusively on " outputs " and Sj exclusively on " inputs ". Sometimes, the 
following notation may be more convenient: 5i x dj = ^iCj. Wc call these operators 
(5e-operators and wc denote their set by A x A, where by A wc denote the set 
of simple 5^ operators. By convention, we define 6e{y{t)} = {1} for each t E Z. 
This means that x de[y{t),u{t)] = Siy{t), dg x 6j[y{t),u{t)] = eju{t). Practically 
speaking, the symbol e corresponds to the empty multiindex {}. Obviously, the 
operators 6^6^ — Si or Se^i = ei, i E Z, are none other than the well-known simple 
shift operators, i.e., SiEeyit) = y(t — i) or Se^iu(t) = u{t — i). The operators 5q, eo, 
give zero delays, i.e. 5Qeey{t) = y(t), 6eeQu{t) = u{t). The 5e-operators of the form 
SoSi^ ■ ■ ■ Si^eoSj^ ■ ■ ■ ej„ are many times called zero terms. 

Example 2.1 Leti= (0, 1,1,2), j = (1, 1, 3, 3, 3), r = (0,0,1,2) ands= (0,1,1,2), 
then diej[y{t),u{t)] = 5(o,i,i,2) e(i,i,3,3,3)[?/(t), w(t)] = 6oSi6i52 eiei€3e3€3[y{t),u{t)] 
= 5oSlS2 ejel[y{t),u{t)] = y{t)y^{t - l)y{t - 2)u'^{t - l)u^{t - 3), furthermore 
Sree[y{t),u{t)] = Sry{t) = 6l5i52y{t) = y\t)y{t - l)y{t - 2) and Sees[y{t),u{t)] = 
esu{t) — eoele2u{t) — u{t)u^{t — l)u{t — 2) 

The order among multiindices implies an order among 5e-operators in a natural way. 
Indeed, we say that ^lej -< Si^ey if either i < i' or i = i' and j < j'. We equip the set of 
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the 5e-operators with two internal operations: the dot-product and the star-product. 
The dot-product corresponds to the usual product among sequences, while the star- 
product corresponds to the substitution of one sequence by another. Specifically, if 
z{t) = 5i£j[y{t),u{t)], s{t) = Syey[y{t),u{t)], i = (ii,i2, • • • ,4), j = (ii,i2, • • • ,Jm), 
i' = (i;,i'2,. . .,4,) J' = ■ • •iJm')- then 6i6j ■6i,ey[y{t),u{t)] = y{t-ii)-y{t-i2) 

■ ■■y{t-in)-y{t-i'i) ■ ■■y{t-i'^,) u{t-ji)-u{t-j2) ■u{t-jjj-u{t~j[) ■ ■■u{t-j'^,). 
Let us now suppose that z{t) = diej[y{t),u{t)], and y{t) = di'w{t), u{t) = eyv{t), 
then Sie^ * [5i', ey][w{t),v{t)] — Sie^[5iiw{t), eyv{t)]. The meaning of the last relation 
is the following: Let us define the maps Z : F x F ^ F,Z[y(t),u{t)] = z(t) = 
6iei[yit),u{t)], Y : F-^ F,Y[wit)] = yit) = dM^), U : F ^ F,U[v{t)] = u{t) = 
ejv{t). The star-product is nothing else than the composition Z o [y, U] = Z[Y, U] 
or by using the other notation diCj * [6i', ey] = Sie^ o ej']= 5iej[5i', ey]. 
The following properties are valid for all the above operations. Their proofs are 
obvious from the definitions and are therefore omitted. 

Proposition 2.1 (a) diSj ■ Si>ey = (5iei'ej©j' (b) S-.e^ * [Si',ey] = S-.r^i^Si^^^^ ■ • -^i'+i^ 
ey+j,ey^j^ ■ ■ ■ ey^j^ (c) Sie^ * {[yi{t) , Ui{t)] + [y2{t) , U2{t)]) ^ * \yx{t),Ux{t)\^8i^ * 

\y2{t)M^)\ 

Part (c) of the above proposition, indicates that the 5e-operators do not behave like 
members of an algebraic ring. 

Example 2.2 Let i = (0,1), j = (1,1,2), i' = (0,0,1) anc? j' = (0,0). Then, 
he; ■8,ey\y{t),u{t)\ = 6odiele2 6l6,el[y{t),u{t)] = d'o6lelele2[yit),u{t)] = y%t)yit - 
l)u'^{t)u^{t - l)u{t - 2). Since i©i' = (0,0,0,1,1) and j©j' = (0,0,1,1,2) we can 

easily verify formula (a) of proposition (2.1). Moreover, Sie^ * [S[r . ey][y{t),u(t)] = 
6oSiele2 * [y\t)y{t - l),u\t)] = y\t)y{t - l)y\t - l)y{t - 2)u\t - l)u\t - 2). Part 
(h), of proposition (2.1), can be verified too, through straight substitution. 

2.3 The ^e-polynomials 

The 5e-polynomials are straightforward extensions of the 5e-operators that are nec- 
essary for the description of discrete non-linear systems with products among input 
and output sequences. These have been introduced and studied in [9], [17]. We 
will go over the main points, adding some new results that will clarify the ideas 
in the present paper. Expressions of the form: A = Z!(ij)eixj Cij<^iej, are called Se- 
polynomials, where by I, J we denote sets of multiindices and Cij real numbers. By 
putting I = {6e}, or J = {eg}, we can transform the above polynomials to pure S 
or e-polynomials. That is A — Z)iei ^ie^i or A = Z^jejCejCj. These are polynomials 
that deal only with input or output signals and produce non-linear polynomial ex- 
pressions of these, without cross-products. By SiA if A is a 5-polynomial, or e^A 
if A is an e-polynomial, we denote the following operations: SiA — Ai^ ■ Ai^ ■ ■ ■ Ai^ 
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with A, = Eiei Cie(^i+,, and e-^A = Aj^ ■ Aj^ ■ ■ ■ Aj„ with Aj^ = EjeJ CejCj+j,. In 
other words we shift all the delays of Ahj ii, then by i2 and so on and finally we 
multiply the obtained polynomials. We will call expressions of the form Ei=i Qe^i 
or J2j=iCej^j, linear 5 or e-polynomials, respectively. The term, which accordingly 
to the order defined previously is ordered highly among the terms of A, is called the 
maximum term of A. By d{A) wc denote the minimum delay of A, in other words 
d{A) = min{min(p),p G I U J}. The largest of the numbers deg{i + j), is called 
the degree of A, denoted also by deg{A). Two 5e-polynomials A and A' are equal if 
A[y{t),u{t)] = A'[y{t),u{t)] for each y{t),u{t) G F and t>0. If A is a (5e-polynomial 
and B, C are a 5 and an e-polynomial, then their star-product, A* [S, C], is defined 
as the composition A o [B, C]. This means that B is replaced into the 5-part of A 
and C into the e-part. Formulae for the calculation of the star-product and addi- 
tional properties, can be found in [9], [10], [17]. These papers prove that the set of 
the pure 5 or 5e-polynomials does not form an algebraic ring. This property makes 
them different to all other similar approaches [1], [3]. The following formulas will 
be used in the proofs of the main theorems of this paper, below. Their proofs are 
straightforward and are therefore omitted. 

Proposition 2.2 (a) Let A he a S-polynomial and B an e-one, then Sie-^ * [A, B] — 

(M)-(eji?) 

(b) Let (JiCj he a d-operator, with i = (ii, ^2, • • • , v)? J ~ {ji^ j^, ■ ■ ■ , jx) and L = 

J2ieii''^i^i> ^ — J2jeJi''TT-j^j linear 5 and e-polynomials, /i, Ji C Z finite sets of 
multiindices, then 

5iej *[L,M]= Yl ^siTs^ ■ ■ ■ rs^rrih^mh^ ■ ■ ■ m/^^^s+ieh+j 
se/f he 

where s = (si, S2, . . . , s^^), h = (/ii, /i2, . . . , h\) and If, the cartesian products 
If = Ii X Ii X ■ ■ ■ X Ii, = Ji X Ji X ■ ■ ■ X Ji. 

V ' ^ . ' 

tp— times X— times 

Property (a) is a remarkable one. It creates, by suitable choosing of n and m, 
any polynomial term including products of certain delays among input and output 
signals. Therefore, it indicates that the star-product acts likewise the derivative in 
the case of continuous systems, (it creates terms with derivatives of any order), [11]. 

Example 2.3 (1) Ji^aei * [25o + ^i, 2eo - ei] = [6162 * {25o + 61)] ■ [el * (2eo - ei)] = 
(251 + ^2) -(252 + (^3) •(2e2-e3)-(2e2-e3) = 16ei5i52 - 16e2e35i52+ Ael5r52 + 9>el5l- 

8e2f,3(^2+ 2e3^2 + 8e2e35i53 + 2el5i5^+ Ael525^ - 46263(52^3 + ^152^:^ 

(2) {25it2 + 5leie2) *[5i -2^2, egei] = 25ie2* [5i -2^2, e^d] +5^6162* [^i -2^2, 60^1] 
= [25-, * {5, - 282)] ■ [62 * e^ei] +[5^ * {5, - 282)] ■ [eie2 * ^e,] = 2(^2 - 25:,) ■ {ele^)+ 
{5i - 252f ■ tlt2 ■ tlt^ = elele3dl+ 2eie3(52+ Aejele^dj- 4ele3d3-4did2elele3 
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The so-called " homogeneous " polynomials are special cases of 5e-polynomials with 
specific applications. 

Definition 2.1 Let 9 = (ii, ^2, • • • , in) be a given multiindex. We call a de-polynomial 
T, "homogeneous" with respect to the multiindex 9 and the constants A,//, if it has 
the form: T = jjI'Seei^ei^ ■ ■ ■ 6^^+ A//"~^5i^ ej^ejg ■ ■ ■ ej„ +A/x""^5i2 ^^h^^ia H 

^V""^ ^i„-i^i„eiiei2 • hA"-V(5n^i2 • --K-i^in H hA"-V(5i2^i3 • ■■K^h + 

A"^ii^j2 • • - (^j^ee = 

= E E AV"-^5i6j (2) 

v=o (ij),iei^,j=6iei 

where X, /j, & R and 1^ is the set of the (p-subindices of the multiindex i, Iq = {e}. 
The following property is of special interest. 

Proposition 2.3 Let A be a 6e-polynomial. If A can be written in the form: 
A — A * [L, M], where L, M are linear 5 and e-polynomials respectively and A is a 
homogeneous Se-polynomial, with respect to an multiindex 9 and the constants A, fx, 
then A^S0*{XL + jiM). 

Proof: Using the 5-symbol, for the delays of L and the e-symbol, for the delays of 
M, we get the following succession: 5g * {XL + ^M) = 5i^5i^ ' ' ' ^i^* (-^-^ + 1^^) = 
5i, * {XL + fiM)- 6i, * {XL + fiM) ■■■ 6,,^* {XL + fiM) = X''{6i, * L) ■ {5i^ * L) 
■ ■ ■ * L) + • • • + X^^I'-\5,, * L) ■ (q, * M) ■ ■ ■ (e,„ * M) + ■ ■ ■ * M) ■ (e,, * M) ■ ■ • 

(e,„ * M) = • • • = X^5i,5i, ■ ■ ■ Si^e, * [L, M] + • • • + Xfi'^-%,e,, ■ ■ ■ e,„ * [L, M] + 

h ii^-Seei^ei^ ■■■^in* [L, M] = [X"-6i^5i^ ■ ■ ■ Si^Se H h A/x^-M.^e,^ • ■ • ei„ H \- 

H^Sgei^Si^ ■ ■ ■ ej„] * [L, M] = A * [L, M], where A is a homogeneous polynomial with 
respect to the multiindex 9 and the constants A, fx. The proposition has been proved. 
□ 

2.4 The Formal JeL- Factorization. 

In this section we shall present the main tool that is used throughout this paper. 
Let {wij} and {sij}. i = 1,2, . . ., j = 1,2, . . . be two sets of parameters, i.e., sets 
of variables that can take arbitrary real values. Let A be a de-polynomial. We 
call 'Formal-5e-Linear-Like- Factorization' of A (F5eL- Factorization for short), any 
expression of the form: 

9 

A = SoSe * [Ag ^, Se] + Se€o * [Se, Ae,i] + E Ck{wij, %)<5i^ejfc * [Lk, Mk] + R 
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where Ck{wij, Sij) are functions of the parameters Wij, Sij, ^, A^^i the 5 and e-hnear 
terms of A, ^i^ej^. G A x A, with min(ifc) = 0, min(jfe) = and 

Lk = WkQ^Q + Wki5i H \-5k , Mfc = Sfeoeo + SkiSi H h 

That is, they are Unear 5 and e-polynomials, where their coefficients are not con- 
stant numbers but the parameters Wij,Sij correspondingly (by the exception of the 
highest delay term). The 5e-polynomial R, with parametrical coefficients, which 
is called the remainder, must contain only zero terms, i.e., terms of the form 
5o5i^ ■ ■ - Si^ eoeji . . . ej^ . The F5eL- Factorization of a given polynomial A is de- 
noted by Formal[A]. For pure ^-polynomials or e-polynomials. wc can restrict the 
above terminology by either ej^ = eg, Mk = eg or Sj^, = 6e, Lk = Se- For instance, 
we say that a (^-polynomial B, has a Formal - (5-Linear - Like - Factorization or, 
simply, an F5L-Factorization, ii B — J2l=iCk{wij)6i^ee * [Lfc,ee] -|- i? or in short 
B = ELi Cfcf^ifc *Lk + R. (See [8] for further details.) 

Example 2.4 A formal 5eL-Factorization of the 5e-polynomial B — 5l + 5iei is 
Formal[B] = 5l*{wio5o +Si) +5oeo*[w2oSo S2oeo + ei]- S2oSoeo*[w3oSo + 5i,eo\ 
-w^2o'^oeo*[^o,-54oeo + ei]+ [-2wio5o5i-wfQ5l+ {W20S20+W30S20 +W2oS4o)Soeo\, where 
Wio, W20, S20, Wso, S40 are the parameters. 

Theorem 2.1 For a given 5e-polynomial A, the F 5 eL- Factorization is unique, as- 
suming that the parameters {wij,Sij} are considered as constants. 

Proof: Wc explain first that by the expression " arc considered as constants " 
we mean that wc treat the parameters like to be specific number. Now, if A has 
only zero terms then the proof of the theorem is trivial and = and R = A. 
Let us suppose that A contains non-zero terms and (^a^ia^ia is its maximum term. 
Using proposition (2.2), (b), we can see that this term also appears in the product 

Ca(Wij,Sij)Si^ej^ * [La, Ma] = Ca{Wij , Sij)Si^ej^* [woaSo + WiaSi + ■■■ + S„,^, Soa^O + 

Sia^i + ■ ■ ■ + Cnal- equating their coefficients we calculate the quantity Ca{wij, Sij) 
uniquely, in effect Ca{wij, Sij) = uja- If we repeat this procedure for the next higher 
order term, we find an expression for the "next" coefficient Ca-i{wij, Sjj). Since this 
expression is a function of Ca{wij, Sij) and of some parameters Wij, Sij considered as 
constants, we conclude that Ca-i{wij, Sij) is also defined uniquely. By induction, we 
can then see that all the coefficients are uniquely determined. The polynomial R 
consists only of zero terms. These terms arise either from the polynomial A or from 
the products Ck{wij, Si^O^i^ej^ * [L^, M^]. The unique determination of the coefficients 
Cjfc entails the uniqueness of R and, thus, the theorem has been proved. □ 

2.5 The Algorithms. 

In what follows we provide the basic algorithmic tool that we will then use to obtain 
a Formal 5eL-Factorization. This is along the lines of Ritt's remainder algorithm 
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[18]. It is based on a kind of division with respect to the star-product [11]. 
The F(^gL-Subroutine. 

Input: A non-linear 5e-polynomial A, without linear or constant terms. 
Initial Condition: The index A = 0. 

DO 

STEP 1: We set A = A + 1. 

STEP 2: Let S = cx{wij, Sij)5iej = cx{wij, Sij)6i^6i^ ■ ■ ■ Si^ej^ej^ ■ ■ ■ ej^ be the 
maximum non-zero term of A, < ii < i2 < ■ ■ ■ < in, < ji < j2 < • • ■ < jm 
and c\{wij,Sij) its coefficient. (Actually, in the first iteration of the algorithm, 
cx{wij,Sij) is always a real number. It then becomes a function of the unknown 
parameters Wij,Sij.) 

STEP 3: We form the linear (5-polynomial 

Lx = wxoSo + wxiSi + WX2S2 H \-Si^ 

where wxo,wxi, ■ ■ ■ are unknown parameters, taking values in R. 
STEP 4: We form the linear e-polynomial 

Mx = sxoeo + sxiei + 8x2^2 H h 

where sxq-, sxi-, ... are unknown parameters, taking values in R. 
STEP 5: We calculate the quantity: 

Rx = ^- cxiwij, Sij)5oSi2-ii ■ ■ ■ (5i„_jieoej2_jj • • • ej^-j^ * [La, Mx] = 

= A- cxiwij, Sij)Si^€^^ * [Lx, Mx] 
STEP 6: We replace the polynomial A with the polynomial Rx- 

UNTIL All the terms of A become zero terms. 
STEP 7: We rename the last value of A as R. 
Output: The quantities Ck{wij, Sy), Si^ej^,Lk, M^, R, k = 1, . . . ,ip. 

The i-VkZ- Algorithm. 
Input: A (5e-polynomial A. 

STEP 1: We decompose A as follows: A = A§ j^+A^^i+Ani, where Ag ^ A^^i, Ani 
are the (5-linear, e-linear and the non-linear parts of A, respectively. 

STEP 2: IF A^i = THEN we set = 0, i? = ELSE using the F6eL- 
Subroutine we take the quantities Ck{wij,Sij),S-n,e^^.,Lk,Mk,R, k = 1, . . . ,ip. 

STEP 3: We set c^+i{wij , Sy) = 1, 6i^^^ej^^^ = SoSe, L^+i = A^ ^, M^+i = 
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Output: The quantities Cfc(wij, Sjj), Ji^ej^, Lyt, Mfc, i?, k = I, . . . ,ip + 2. 

Theorem 2.2 The FSeL- Algorithm terminates after a finite number of steps. If 

Ck, Si^^ej^, Lk, Mk, R, k = 1, . . . , ip + 2, are its outputs, then the quantity Formal[A] — 
Ylk=i (^-k^h^jk * [LkjMk] + R is the Formal - SeL - Factorization of the given de- 
polynomial A. 

Proof: If we follow the above algorithm step by step, we get Ri = A—ci{wij, )5i^ejj* 
[Li,Mi]. By means of proposition (2.2), (b), we conclude that this operation elim- 
inates at least one non-zero term of A. Next, we get R2 = Ri — C2{wij, Sij)di,^dj2 * 
[L2,M2], which eliminates another non-zero term of A and so on. These repeated 
operations, ensure that all the non-zero terms of A will be finally eliminated and 
therefore the procedure will terminate. This also means that the remainder R will 
contain only zero terms. Furthermore, for the operators ^i^^ej;^ in Step 5 of the 
F(5eL-Subroutine, we have min(iA) = min(jA) = 0. Additionally, by reverse substi- 
tution we get A = Ci{wij, s.y)5ijej^ * [Li, Mi] + C2{wij, Sij)6i^6j^ 
— YX=i ^ki.'^ij) ^ij)^ik^jk * [-^k, Mk] + R. Finally, all the above ensure that the out- 
put polynomials Si^ej^, L^, M^, R, k — 1, . . . ,ip + 2 and the coefficients Ck{wij, s^), 
constitute a formal (5eL-Factorization for the given 5e-polynomial A. □ 



2.6 Evaluations 

Let A be a 5e-polynomial and Formal[A], its Formal (5eL- Factorization. If the 
parameters Wij,Sij take values according to a set of substitution rules U, we say 
that the Formal[A] is evaluated over the set U, thus writing: Formal[A] ^. More 
rigorously, let W = (wij) be the set of the u'-parameters, written as a vector and 
S = (ski) the vector of the s-parameters. r = (ajj), q = {bki) are vectors of real 
numbers, which are in one-to-one correspondence with the vectors W, S. We say 
that these parameters follow the rule (r, q), thus writing (W, S) (r,q), if the 
following substitutions are valid {wij, s^i) = (ajj, bki). Let N, Q be two sets of rules, 
N = {ri, r2, . . . , ta}, Q = {qi, q2, . . . , q^} and U C N x Q, then 

« r 3 

Formal[A] = IJ 1 2]] Ck{wij, Shi)6i^ej^ * [Lk, Mk] + R, 

a=l U=l 

wtth, (W, S) ^ (r, q)^ G U} 

The set of substitution rules, U, may be finite or infinite. Often we denote this by 
U = {{wij,Ski) = {aij,bki)a,cr = 1,2,...}, or U = {{wij = aij,Ski = bki), {wij = 
KjJ ^kl = b'ki), •••,}• 
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Example 2.5 Let B and Formal[B\ he as in example (2.4)- The vectors of the pa- 
rameters are W = {wio, W20, W30), S = (s2o, S40). For instance, let us take the follow- 
ing sets of rules N = {(1, -1, 2)}, Q = {(-1, i)}, then U = {(wio, W20, W30, S20, S40) 
= (1, -1, 2, -1, 1)} anc? Formal[B] ^ = dl* {Sq + Si) +5Qeo*[-So + 5i, -eo + ei]- 

{-l)5oeo*[2So + 5i,eo\ -(-l)5oeo * [(^o, |eo + ei] +{-2SoSi - 5l - ISoeo). The rules 

N = {{0,ip,-Lp{e + 1))}, Q = {(1,^)}, will give U = {(u^io, w^2o, «^30, S20, S4o) = 
(0, (fi, —(p{6+l), 1, 9), (/?, 6* e R and Formal[B] ^ — 5l*di -\-doeo*[ip5o-\-Si^ ^o+^i] — 



5oeo* [-(p{9+ l)So + Si,eo] -ipSoeo* [So,9eo + ei]. 

3 Part II-The Simplification-Linearization Algo- 
rithms 

In this section wc present the simphfication procedures of our methodology. Before 
we examine them, we would like to begin by a short algebraic description of the 
non-linear discrete input-output systems, via the notion of 5e-polynomials. 

3.1 The (5e-polynomials and Non-linear Discrete Systems 

Suppose that we have a system of the form (1). By using de-operators we can rewrite 
this as follows: 

(^Ayi't) + 13 13 aij6i6jy{t) + ^ I] • • • 13 (^iii2...iJiA^ ■ ■ ■ Si„y{t) = 

= ^ keiuit) + hjeiejuit) + ■ ■ ■ + ^ ■ ■ ■ ^ ■ ■ ■ ^i^e^, ■ ■ ■ €j^[y{t),u{t)] (3) 

We can describe the above relation shortly, by writing Ay{t) = Bu{i) + C[y{t),u{t)], 
where A is a 5-polynomial, B an e-polynomial and C a pure 5e - one. The causality 
and solvability of the system is guaranteed by the inequality: d(Ai) — d{A) < 
mm{d{B),d{C)}, where Ai is the linear part of A. This inequality means that the 
lower delayed term of the output (i.e. y{t)) appears in the linear part of the system 
and thus we can solve (1) with respect to y{t) in a direct way. 

To each nonlinear system of the form (1) we assign a vector of real numbers 
Yo — {Vo, Ui,---: Vk-i) which gives the so called initial conditions: y{0) — y^, y{l) — 
yi, . . . , y{k — 1) = yk-i, where k is the maximum delay appeared in the output 
signal. Since the signals, involved in (1), are causal, i.e. y{t) = 0, u{t) = 0, for 
t < 0, for each given vector of initial conditions any input signal u{t) determines a 
unique output signal that satisfies (1), for n > k. Let us have two causal nonlinear 
systems: Ay{t) = Bu{t) + C[y{t),u{t)\ with yo = (yo, yi, ■ ■ ■ , 2/fc-i) and A'^{t) = 
Bv{t) + C['ip{t), v{ty\ with ifjQ = {ipo, -01, ... , i^x-i) and k > X. We say that the two 
systems " operate " under the same initial conditions or that their initial conditions 
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are identical if = i/o, ipi = iji, . . ., VA-i = Ux-i and yx = ip{X), yx+i = ipiX + 1), 
. . ., yk-i = tp{k — 1). In other words we must give as initial conditions to the system 
which starts to product outputs later, the corresponding outputs of the other system. 
Two systems Aiy{t) = B^u{t) + Cr[y{t),u{t)] and A2y*{t) = B2v{t) + C2[y*{t),v{t)] 
are equivalent, if y{t) = y*{t), whenever u{t) = v{t) and initial conditions are 
identical. A " simplification " method for a given non-linear system consists in 
discovering systems that are equivalent to the original one but have a '' less '' complex 
structure. This procedure allows us to replace, if necessary, the original system with 
the simpler one. 

A 5e-polynomial is called 'proper' if the following two facts obtain: i) The minimum 
delay of A, a = d{A), appears only in the 5-part of certain terms of A. ii) The power 
of 6a, in the above terms, is equal to one. Obviously, all polynomials with linear 
parts that contain the lowest delay terms are proper. The following theorem is an 
extension of a useful result that appears in [13]. The proof is along similar lines. 

Theorem 3.1 The equality A[y{t),u{t)] = A[y*{t),u{t)], \/t e Z, implies that 
y{t) — y*{t), provided that A is proper and y{t),y*{t) have the same initial con- 
ditions. 

A direct consequence of the above theorem is the following result: 

Corollary 3.1 Let A be a proper 5-polynomial, then the equality Ay[t) = 0, Vt e Z 
implies that y{t) = under zero initial conditions. 

All the above terminology can be applied to systems without cross-products, in 
a direct way. These are systems of the form: 

aiy(t-i)+J2 Yl (^ijyit-'>')y(t-j)^ — ^Yl ■■■Yl aiii2...iny(t-^i)y(t-i2) ■ ■ ■ y(t-in) = 

= Yl - + Z] D biju{t - i)u{t - j) H h X] • • • Z] ""(^ - ii) • • • "(^ - jk) 

Using 8 and e-operators we rewrite the above system as 

aidi + YY. (^iA^j + ■ ■ ■ + H ■ ■ ■ H ai^h-in^iAi ■ ■ ■ Sin) y{t) = 
= {Y bi^i + YY1 K^i^j + • • • + • • • S • • • 

or shortly Ay{t) = Bu{t), where A a 5-polynomial and B an e-polynomial. The 
causality is guaranteed by the inequality: d{Ai) — d{A) < d{B), where Ai is the 
linear part of A. Two systems Aiy{t) = Biu{t) and A2y*{t) = B2v{t) arc equivalent, 
if y{t) = y*{t), whenever u{t) = v(t) and initial conditions are identical. A " 
simplification " procedure for these systems is defined by the same way, as before. 
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3.2 The Case without Cross-Products 



In this section we are dealing with systems which do not contain products among 
input and output signals. Actually, we want to simplify them via the use of a 
symbolic algorithm. Clearly, of all such 'simpler' systems, linear systems are the 
most desirable ones. Accordingly, in this section we ask under what circumstances a 
non-linear system of the form Ay{t) — Bu{t) is equivalent to a linear one. In other 
words, wc arc looking for a linear system Aiy*'{t) = Biv{t), such that y*{t) = y{t), 
whenever u{t) = v(t). We are not involved with stability questions. The algorithm 
upon discussion is: 



The Linear-Equivalence Algorithm I 



Input: The nonlinear d, e-polynomials A and B. 

STEP 1: By means of the F(5eL-Algorithm, we find the quantities: 

UJ 

Formal[A] = ^ Ck{wij)6i^^ *Lk + R§ 
k=l 

f 

Formal[B] = ^ c^(sij)ej^ *Mp + Re 
p=i 

STEP 2: We form the set of rules: U = {{wij,Sij) = {aij,bij),aij,bij G R}, 
such that the following are valid simultaneously 



= Re 

u 



= 

u 



• The sets L = {L : L = gcd( Lfc )} and M = {M : M = gcd( Mp )} 
are not void, i.e. L / and M / 0. 

• The polynomials A, B, where yl = 1 * L, L G L, 5 = 5 * M, M G M, are 
proper and A = B. 

STEP 3: IF [/ / THEN goto the output ELSE the method fails. 



Output: The quantities L,M 

Theorem 3.2 Let us suppose that we have the causal nonlinear discrete system 
Ay{t) — Bu{t), A, B 5 and e - polynomials with d{A) < d{B). Ifh and M, are 

the outputs of the previous algorithm, then each linear system of the form Ly{t) — 
Mu{t), L e L, M e M is equivalent to the original nonlinear system. 
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Proof: Let, W and S, be the sets of parameters, appearing in the formal factor- 
ization of A, and B, during the run of the algorithm. The existence of the solution 
implies that C/ 7^ and therefore we can find at least, one pair of vectors, r = (%), 
q = such that the rules W — r, S ^ q make all the relations of the step 2, 
true. Hereafter, we shall consider that all the polynomials, we work with in our 
proof, have been evaluated coherently to the above rules. We have now, for the 
polynomials A and B: 

^=Y1 ^k{wij)6i^ * Lk + R^ = Y1 Ck{wij)Sii^ *Lk*L + Rg^A*L + R^ 

k=l k=l 



B=J2 Cpi^ijhp *Mp + Re = Yl Cp(sii)ej, * * M + Re ^ B * M + Re 
p=i p=i 

Taking into account the fact that the equalities of the step 2 are true, we can 
easily get that A — B, R^ — Re — 0. The original system can now be re- written 

as i * Ly{t) = B * Mu{t). We define yi{t) = Ly{t), Ui{t) = Mu{t) and thus 
Ayi{t) = Bui{t). Since A = B and A,B proper, we conclude that yi{t) = Ui{t), 
which means that Ly{t) = Mu{t). To finish the proof we have to show that this 
linear system is causal too. In other words, we must prove that d{L) < d(J\d), but 
from the causality of the original system we get d{A) < d{B) and d{A) — d{A)-\-d{L), 
d{B) = d{B) + d{M). A = 13 implies the desired result. All the above statements 
hold on in the case of any other pair of rules (p,y), belonging to U and thus, any 
set of polynomials L G L, M G M can form an exact linear-output linearization of 
the system Ay{t) = Bu{t). The proof has been completed. □ 

Remark 3.1 The main advantage of the above procedure is that it gives a set of 
solutions and it works when the polynomials A, B have no linear parts. Nevertheless, 
if the polynomials A and B do have linear parts, then two facts can be easily proved: 
(a) the properness condition of step 2 is not needed (b) the polynomials L and M, 
are factors of the linear parts of A and B. 

Example 3.1 To clarify our ideas we present the following examples: 
(a) Let us consider the nonlinear discrete system: 

Ay{t)y{t - 1) + 2y{t)y{t - 2) + 22y\t - 1) + 21y{t - l)y{t - 2) + hy\t - 2) = 

= u{t - l)u{t - 2) + 2u^{t - 2) - 3M(t - l)u{t - 3) - 2lu{t - 2)u{t - 3) + 45«2(t - 3) 
or using 6 — operators: 

(ASoSi + 28062 + 2251 + 2151^2 + ^5l)y{t) = 
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Repeats 


Non-linear Output 


Linear Output 


5 


0.158373 


0.158373 


100 


-0.34 6706 


-0.34 6701 


300 


0.024609 


0.024609 



Table 1: Simulation of the Simple Linear Simplification. 

= (eiea + 2el - Seieg - 21e2e3 + ^^el)u{t) 

shortly Ay{t) = Bu{t). Applying the Linear-Equivalence Algorithm I, we take the 
following factorizations: 

A = 5^0 * ('"^10(^0 + Wix5i + ^2) + (21 - 1{)wu)5q5i * (^20^0 + 

+ (22 - bwl^ - 21«;2o + muiiW2o)5l * (w^o^o + ^^i) + 

B = 45eo * (sioeo + snei + Si2e2 + €3) + (-21 - 90si2)eoei * (s2oeo + S2iei + €2)+ 

+(-3 - 90sii + 2IS21 + 90si2S2i)eoe2 * (saoCo + ei) + Cg^oei * (ssoCo + 

+(2 - 45si2 + 2IS21 + 90si2S2i)eo * (s4oeo + S4iei + 62) + d^el * (seoCo + ^i) + Re 

(The terms Cg, Cg, are not presented explicitly, due to their huge size.) The set 
U , that is the set of the values of the parameters Wij, Sij, which satisfy the relations 
of Step 2, is the following: U = {{wio = 0,wii = 2,W2o = 2,t«3o = 2, Sio = 
0, sii = 0, S12 = -|, S20 = 0, S21 = -|, S30 = 0, S40 = 0, S41 = -|, S50 = 0, S60 = 0), 

{Wio = 0,Wu = 2,W20 = 2,W30 = 2, Sio = 0,Sii = 0,Si2 = -|,S20 = 0, S21 = 

S30 = — |, S40 = 0, S41 = — |, S50 = — |, Sqo = — |)}- These values will give A — B — 
55^ + SqSi, which are proper polynomials, and L — 25o + Si, M — €1 — 3e2. Hence 
[/ 7^ and the sets L, M are L = {25o + ^i}; M — {ei — 3e2}- This means that the " 
linearization " of the above system is the linear system: {2SQ-\-6i)y{t) — (ei — 3e2)w(t) 
or 

2y{t) + y{t - 1) = u{t - 1) - ?>u{t - 2) 

To verify that the two systems, the nonlinear and the linear one, have the same dy- 
namic behavior we present, in table 1, their responses under the same input sequence 
u{t) — rnd{l), and the same initial conditions, 
(b) We have the causal nonlinear discrete system: 

2y{t) + 2y{t - 1) + ^y{t - 2) + y{t - l)y{t - 2) - ^y\t - 2) - ^y(t - l)y{t - 3)- 
-^Z/(^-2)y(t-3) + ^|/2(t-3)-^2/(t-l)|/(t-4) + ^y(t-2)|/(t-4) + i|/(t-3)t/(i-4) = 
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= 2u{t - 1) + u{t - 2) +u{t - 2)u{t - 3) - u^{t - 3) - u{t - 2)u{t - 4) + u{t - 3)n(t - 4) 

Using the same methodology, we can find that the linear system: {5q + \8i)y{t) — 
5iu{t) or 

y{t) + \y{t-i) = u{t-i) 

consists the simplification upon request. We observe here that the polynomials A 
and B have linear parts Ai = 2Sq + 26 1 + 1^2, Bi = 2ei + €2 and therefore the linear 
polynomials L = Sq + ^5i, M = 61 are factors of Ai and Bi (Remark 3.1). 



3.3 The Cross-Products Case 

Let us suppose now that we have a causal non-hnear discrete system with cross- 
products of the form Ay{t) = Bu{t) + C[y(t), ^(t)]. As we mentioned before, a 
hnearization method for this system consists in discovering hnear systems that are 
equivalent to the original one. The algorithm for this category of systems is: 



The Linear-Equivalence Algorithm II 

Input: The (5-polynomial A, the e-polynomial B and the (5e-polynomial C. 

STEP 0: We define the (5e-polynomial As^s A = A - B - C.{\Ne transfer all 
the terms at the left-hand side of (3)). 

STEP 1: By nneans of the F(5eL-Algorithm, we find the Formal (5eL-Factorization 
of the polynomial A: 

a 

A=Y, (^kiwij, Sij)6i^ej^ * [Lk, Mk] + Ra 
k=l 

STEP 2: We define the set of rules U — {{wij,Sij) — (aij,hij), ajj,6yR} such 
that the following are valid simultaneously 

• Rj^ u~^' eliminate the remainder). 

• The sets L = {X : L = gcd ( Lk |^) }, M = {M : M = gcd ( Mk ^)}, 
k = are not void. (We find the common factors of the inear 
polynomials). 

• The (5e-polynomials, Ap, p = \,...,k, where A = (X]p=i «piip) * [L,M], 
L h, M £ M, Op G R, are homogeneous, with respect to the multiindices 
ip,p = l,...,k and the standard constants A,//. (We factorize A with 
respect to the common factors). 

• The polynomial A = Ylp=i0.pSip< is proper. (The multiindices i^ are the 
multiindices of the homogeneous polynomial). 
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STEP 3: IF [/ ^ THEN goto the output ELSE the method fails. 



Output: The quantities L,M, A, 
The condition U ^ ^ is equivalent to the solvabihty of an algebraic system of 
polynomial equations. Indeed, since the remainder Rj^ is a 5e-polynomials with 

coefficients which involve the parameters {wij,Sij), the condition Rj[ ^ = 0, is 
equivalent to a system of polynomial equations. This system can be solved by 
successive substitutions. The first equation contains only one unknown parameter 
and thus it is solvable. Substituting the value will obtain to the second equation 
we take a solution for the second parameter and so on. Among the solutions of this 
algebraic system, if any, we can choose those particular one, which can satisfy and 
other conditions, like minimality in the number of evaluated parameters and so on. 

Theorem 3.3 Let Ay (t) = Bu{t)+C[y{t) , u{t)], A a 6 -polynomial, B an e-polynomial 
and C a de - polynomial, be a causal non-linear discrete system that contains prod- 
ucts among input and outputs. //L,M,A; fx, are the outputs of the previous algo- 
rithm, then, any linear system of the form: XLy{t) — —nMu{t), L e L, M e M is 
equivalent to the original non-linear system. 

Proof: Let W = {wij} and S = {sjj}, be the sets of the parameters that appear in 
the Formal 5eL-Factorization of A, during the run of the algorithm. The existence 
of the solution entails that U ^ ^ and, therefore, that there is at least one pair 
of vectors, r = (a.y), q = (fey), such that the rules W — r, S — q make all the 
relations of the step 2, true. Hereafter, we will assume that all the polynomials in 
the proof, have been evaluated and cohere with the above rules. For the polynomial 
A we now have: 

9 f ^ ^ 

k=l k=l 

k 

= (E «pip) * [L, M] + Ra 
p=i 

By taking into account the fact that the equalities in step 2 are true, we can easily 
see that Rj^ — 0. Moreover the polynomials Ap are homogeneous. This means that 

taking into account proposition 2.3 we get A — Opf^ip) * (XL + jiM), or A = A* 

{\L + jiM). The original system can be now re- written as A* {\L + iiM)[y{t),u(t)] = 
We define uj{t) = XLy{t) + fiMu{t), and the system becomes Auj{t) = 0. The last 
condition of the step 2, guarantee that A is proper. This means that u{t) = 
(corollary 3.1 ) or XLy{t) + iiMu{t) — and thus XLy{t) — —iiMu{t). To conclude 
the proof we need to show that this linear system is also causal. In other words. 
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we must prove that d{XL) < d{—fiM). From the causahty of the original system 
we get: d{Ay) < d{A^), where d{Ay) is the minimum delay of the output signal at 
the polynomial A, d{Ay) is the minimum delay of the input. Additionally, using the 
equalities, presented above, we get: d{Ay) — d{A) + d{L), d{Ay) — d{A) + d{M). 
The elimination of d{A) entails the desired result. All the above statements hold for 
any other pair of rules (p, y) belonging to U and thus, any set of polynomials L Eh, 
M e M can form an exact linear-output linearization of the nonlinear system. This 
completes the proof. □ 

Example 3.2 We have the system: 

6y(t) - ny(t - 1) + 6y(t - 2) - y(t - 3) + 36y(t - l)y(t - 2)- 

-30y{t - l)y{t - 3) - 30y\t - 2) + 25y{t - 2)y{t - 3) - y{t - 3)y{t - 4) = 

= u{t - 1) + u{t - 2) - 2u{t - 3) - u{t - 2)u{t - 3) + 4u{t - 3)u{t - 4)+ 

+6y{t - l)u{t - 3) - 5y{t - 2)u{t - 3) + 6u{t - 2)y{t - 2)- 

-by{t - 3)u{t - 2) - 2u{t - 3)y{t - 4) - 2y{t - 3)u{t - 4) 

Putting all the terms at the left-hand side of the above equation and using 5,e and 
5e-operators we find A. That is: 

A = 66o - lUi + 662 - ^3 - ei - es + 2e3 + SG^i^s - 3O6163- 

-6^163 - 30^2 + 25^2^3 + 552^3 - 6(^262 + 55^62+ 

+€263 - 46364 + 26364 + 28364 - 8364 

Following the Linear- Equivalence II Algorithm step by step we get U ^ for a 
specific set of values for the parameters, denoted as (r, q). (For reasons of brevity, 
we have not shown all calculations in detail). Using the substitution rule (W, S) 
(r, q), we rewrite the polynomial A as follows: 

A = 5o6e * [6S0 - 55i + 82, -ei - 262] + 56q * [65o - 55i + ^2, -ei - 262]- 

—5iee * [6^0 - 5(^1 + 52, -ei - 262] - 5e6i * [6^0 - 5(5i + 82-, -ei - 2e2] + 
+5i52ee * [6^0 - 5^1 + 82, -ei - 262] + 8162 * [Q8q - 55i + ^2, -ei - 2e2] + 
+^261 * [65o - 5(5i + ^2, —61 — 262] + 8eei62 * [65o - 55i + ^2, -ei - 262] = 

= [{8Q6e + ^eCo) " (^l^e + 8e6i) + {8i82ee + 8162 + 5261 + 

+5eeie2)] * [6^0 - 5(5i + 82, -ei - 2e2] 

The polynomials in the parenthesis on the right of the above expression, are homoge- 
neous polynomials in relation to the multiindices, 0, 1, (1,2) and constants X — /i — 
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Repeats 


Non-linear Output 


Linear Output 


50 


0.922726 


0.922726 


500 


0.656207 


0.656207 


1000 


1.15928 


1.15928 



Table 2: Simulation of the Linear Equivalence 

1. We see that = = 1; 0,1 = W = —1; Q'(i,2) = ^(1,2) = 1 o-nd therefore we 
construct the polynomial A = — 5i + 5i52, which is proper. Thus, finally we get 
L = {6^0 — 55i}, M = {— ei}, \ = 11 — 1. Therefore the desired linear system is 
Ly{t) = -Mu{t) or (6(5o - 55i + 52)y{t) = (2e2 + ei)u{t) or 

6y{t) - 5y{t - 1) + y{t - 2) = u{t - 1) + 2u{t - 2) 

Simulations are presented in Table 2. We took u{t) — rnd{l) as input. 

4 Concluding Remark 

The aim of this paper was to describe algebraic computational methods for the 
simplification of a general class of non-linear discrete input - output systems that 
contain products between input and output signals. Actually, we developed an 
approach to the Linear - Equivalence problem. The entire approach is based on 
a proper framework that involves the so-called 5e-operators and the star-product 
operation. We hope to be able to present current work on further applications of 
this method to more concrete questions in a future paper. 
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